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Abstract 

The Higgs sector of the Minimal Supersymmetric Model (MSSM) is a 
CP-conserving two-Higgs doublet model that depends, at tree-level, on two 
Higgs sector parameters. In order to accurately determine the phenomeno- 
logical implications of this model, one must include the effects of radiative 
corrections. The leading contributions to the one-loop radiative corrections 
are exhibited; large logarithms are resummed by the renormalization group 
method. Implications for Higgs phenomenology are briefly discussed. 

1 Introduction 

The Standard Model with minimal Higgs content is not expected to be the 
ultimate theoretical structure responsible for electroweak symmetry breaking [0,0. 
If the Standard Model is embedded in a more fundamental structure characterized 
by a much larger energy scale [e.g., the Planck scale, which must appear in any 
theory of fundamental particles and interactions that includes gravity), the Higgs 
boson would tend to acquire mass of order the largest energy scale due to radiative 
corrections. Only by adjusting {i.e., "fine-tuning") the parameters of the Higgs 
potential "unnaturally" can one arrange a large hierarchy between the Planck 
scale and the scale of electroweak symmetry breaking |]^,^. The Standard Model 
provides no mechanism for this, but supersymmetric theories have the potential to 
address these issues. In a supersymmetric theory, the size of radiative corrections to 
scalar squared-masses is limited by the exact cancelation of quadratically divergent 
contributions from loops of particles and their supersymmetric partners. Since 
supersymmetry is not an exact symmetry at low energies, this cancelation must be 
incomplete, and the size of the radiative corrections to the Higgs mass is controlled 
by the extent of the supersymmetry breaking. The resolution of the naturalness 
and hierarchy problems requires that the scale of supersymmetry breaking should 
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not exceed 0(1 TeV) 0. Such "low-energy" supersymmetric theories are especially 
interesting in that, to date, they provide the only theoretical framework in which 
the problems of naturalness and hierarchy can be resolved while retaining the Higgs 
bosons as truly elementary weakly coupled spin-0 particles. 

The Minimal Supersymmetric extension of the Standard Model (MSSM) con- 
tains the Standard Model particle spectrum and the corresponding supersymmetric 
partners [§|,[7|. In addition, the MSSM must possess two Higgs doublets in order to 
give masses to up and down type fermions in a manner consistent with supersym- 
metry (and to avoid gauge anomalies introduced by the fermionic superpartners 
of the Higgs bosons). In particular, the MSSM Higgs sector is a CP-conserving 
two-Higgs-doublet model, which can be parametrized at tree-level in terms of two 
Higgs sector parameters. This structure arises due to constraints imposed by su- 
persymmetry that determine the Higgs quartic couplings in terms of electroweak 
gauge coupling constants. 

In section 2, I review the general structure of the (nonsupersymmetric) two- 
Higgs-doublet extension of the Standard Model. By imposing the constraints of 
supersymmetry on the quartic terms of the Higgs potential (and the Higgs-fermion 
interaction) one obtains the Higgs sector of the MSSM. The tree-level predictions 
of this model are briefly summarized in section 3. The inclusion of radiative cor- 
rections in the analysis of the MSSM Higgs sector can have profound implications. 
The most dramatic effect of the radiative corrections on the MSSM Higgs sector 
is the modification of the tree-level mass relations of the model. The leading one- 
loop radiative corrections to MSSM Higgs masses are described in section 4. These 
include the full set of one-loop leading logarithmic terms, and the leading third 
generation squark- mixing corrections. In section 5, the leading logarithms are re- 
summed to all orders via the renormalization group technique. A simple analytic 
formula is exhibited which serves as an excellent approximation to the numerically 
integrated renormalization group equations. Numerical examples demonstrate that 
the Higgs masses computed in this approximation lie within 2 GeV of their actual 
values over a very large fraction of the supersymmetric parameter space. Finally, 
some implications of the radiatively-corrected Higgs sector are briefly explored in 
section 6. Certain technical details are relegated to the appendices. 



2 The Two-Higgs Doublet Model 

I begin with a brief review of the general (non-supersymmetric) two-Higgs dou- 
blet extension of the Standard Model 0. Let $i and $2 denote two complex Y = 1, 
SU(2)^ doublet scalar fields. The most general gauge invariant scalar potential is 
given by 
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+ iAi(<l>I$i)2 + \\2{^\^2f + A3($I$l)($|$2) + A4($I$2)($^$l) 

+ {iA5($l$2)'+ [A6($I$l) + A7($|$2)]$l$2 + h.C.} . (1) 

In most discussions of two-Higgs-doublet models, the terms proportional to Ag and 
A7 are absent. This can be achieved by imposing a discrete symmetry $1 — »• — $1 
on the model. Such a symmetry would also require mi2 = unless we allow a soft 
violation of this discrete symmetry by dimension- two terms.Q For the moment, I 
will refrain from setting any of the coefficients in eq. (|l]) to zero. In principle, m^2, 
A5, Ae and A7 can be complex. However, for simplicity, I shall ignore the possibility 
of CP-violating effects in the Higgs sector by choosing all coefficients in eq. (|l|) to 
be real. The scalar fields will develop non-zero vacuum expectation values if the 
mass matrix m'fj has at least one negative eigenvalue. Imposing CP invariance and 
U(l)p.(^ gauge symmetry, the minimum of the potential is 



*i)=-^ . . (*2)=-^ : . (2) 





where the Vi are assumed to be real. It is convenient to introduce the following 
notation: 



v' = vl + vl = ^ = (246 GeV)2 , t = tan/3 = ^ . (3) 

5 Vl 



Of the original eight scalar degrees of freedom, three Goldstone bosons {G^ and 
G") are absorbed ("eaten") by the and Z. The remaining five physical Higgs 
particles are: two CP-even scalars {h^ and H^, with m/jo < niHo), one CP-odd 
scalar (A^) and a charged Higgs pair [H^). The mass parameters mn and 777,22 can 
be eliminated by minimizing the scalar potential. The resulting squared masses 
for the CP-odd and charged Higgs states are 



2 

9 2 



m~jj± = 777^0 



' +yiX,-X,). (4) 



-'^This latter requirement is sufficient to guarantee the absence of Higgs- mediated tree- level 
flavor changing neutral currents. 
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The two CP-even Higgs states mix according to the following squared mass matrix: 

^ ^2 Xicl + 2\spcp + (A3 + U)S0C0 + XqcI + AyS^ \ 

V (A3 + XA)sj3Cp + Aec^ + -^7^ A24 + 2\jSi3Cj3 + A5CI y 

where sp = sin /3 and = cos /3. The physical mass eigenstates are 
= (v^Re$? — vi) cos a + (\/2Re$^ — ^2) sina , 
= -(V2Re$? -ui)sina+ (\/2Re $^ - ^2) cosa . (6) 

The corresponding masses are 

<o,,o = I + Ml, ± V(M?i-A^i2)' + 4(M?2)^ 
and the mixing angle a is obtained from 



(7) 



sin 2a 



2Mj, 



{Ml, - Ml,r + 4{Mhy 



Mil - M22 



cos2q; = . (8) 

^{Mj^-Mhy + ^Mh)' 

The phenomenology of the two-Higgs doublet model depends in detail on the 
various couplings of the Higgs bosons to gauge bosons, Higgs bosons and fermions. 
The Higgs couplings to gauge bosons follow from gauge invariance and are thus 
model independent. For example, the couplings of the two CP-even Higgs bosons 
to W and Z pairs are given in terms of the angles a and /3 by 

9h°vv = dv'^v sin(/3 - a) 

9h0vv = 9v^W cos(/5 - a) , (9) 

where 

There are no tree-level couplings of or H"^ to VV. Gauge invariance also 
determines the strength of the trilinear couplings of one gauge boson to two Higgs 
bosons. For example, 

g cos(/9 — a) 



9h°A°Z — 



2 cos 6w 



2cos^^ ■ ^^^^ 
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In the examples shown above, some of the couphngs can be suppressed if either 
sin(/3 — a) or cos(/5 — a) is very small. Note that all the vector boson-Higgs boson 
couplings cannot vanish simultaneously. From the expressions above, we see that 
the following sum rules must hold separately for V = W and Z: 

2 2 2 2 

9hOVV + 9hOVV = 5'y"^y ; 

2 2 9^ 
9hOAOZ + 9 HO AO Z = ■ (■'-^) 

4 COS'' Ow 

These results are a consequence of the tree- unit arity of the electroweak theory 0. 
Moreover, if we focus on a given CP-even Higgs state, we note that its couplings to 
VV and A'^V cannot be simultaneously suppressed, since eqs. (|9|)-(pJD imply that 



9lzz + 4^^z91aoz = ^ Ta ' ^^3) 
COS"' Ow 

for h = or H^. Similar considerations also hold for the coupling of and 
to W^H"^. We can summarize the above results by noting that the coupling of 
and to vector boson pairs or vector-scalar boson final states is proportional to 
either sin(/3 — a) or cos(/9 — a) as indicated below 0,0. 

cos(/? — a) sin(/9 — a) 

H^W+W- h^W+W- 

H^ZZ h^ZZ 

ZA'^h'' ZA^H^ (14) 

ZW^H^h^ ZW^H^H^ 

Note in particular that all vertices in the theory that contain at least one vec- 
tor boson and exactly one non- minimal Higgs boson state {H^, A^ or are 
proportional to cos(/5 — a). 

The 3-point and 4-point Higgs self-couplings depend on the parameters of the 
two-Higgs-doublet potential [eq. (|I|)]. The Feynman rules for the trilinear Higgs 
vertices are listed in Appendix A. The Feynman rules for the 4-point Higgs vertices 
are rather tedious in the general two-Higgs-doublet model and will not be given 
here. 

The Higgs couplings to fermions are model dependent, although their form is 
often constrained by discrete symmetries that are imposed in order to avoid tree- 
level flavor changing neutral currents mediated by Higgs exchange |rT| . An example 



of a model that respects this constraint is one in which one Higgs doublet (before 
symmetry breaking) couples exclusively to down-type fermions and the other Higgs 
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doublet couples exclusively to up-type fermions. This is the pattern of couplings 
found in the MSSM. The results in this case are as follows. The couplings of the 
neutral Higgs bosons to // relative to the Standard Model value, gnif /2mw, are 
given by (using 3rd family notation) 



— Sill Of 

h^bb : = sinip — a) — tan P cos(p — a) 

cosp 



h'^tt : _ gjj^(^^ _ Q,) _|_ cos(/3 — a) 

sin p 



_ cos oc 

H bb : = cos(/5 — a) + tan /3 sin(/5 — a) 

cosp 



— sill cx 

H^tt : = cos(p — a) — cot P sin(/3 — a) , 

sm p 

A'^bb : 75 tan /? , 

: 75 cot/?, (15) 

(the 75 indicates a pseudoscalar coupling), and the charged Higgs boson coupling 
to fermion pairs (with all particles pointing into the vertex) is given by 

9H-fb = ^ttI [mtCot/9(l + 75)+mfetan/?(l-75)]. (16) 

The pattern of couplings displayed above can be understood in the context of 



the decoupling limit of the two-Higgs-doublet model [|12| , p!3| . First, consider the 
Standard Model Higgs boson (0°). At tree- level, the Higgs self-coupling is related 
to its mass. If A is the quartic Higgs self-interaction strength [see eq. (^7])], then 
A = m^o/f^. This means that one cannot take m^o arbitrarily large without the 
attendant growth in A. That is, the heavy Higgs limit in the Standard Model 
exhibits non-decoupling. In models of a non-minimal Higgs sector, the situation is 
more complex. In some models (with the Standard Model as one example), it is not 
possible to take any Higgs mass much larger than 0{v) without finding at least one 
strong Higgs self-coupling. In other models, one finds that the non-minimal Higgs 
boson masses can be taken large at fixed Higgs self-couplings. Such behavior can 
arise in models that possess one (or more) off-diagonal squared-mass parameters in 
addition to the diagonal scalar squared-masses. In the limit where the off-diagonal 
squared-mass parameters are taken large [keeping the dimensionless Higgs self- 
couplings fixed and ^ C>(1)], the heavy Higgs states decouple, while both light 
and heavy Higgs bosons remain weakly-coupled. In this decoupling limit, exactly 
one neutral CP-even Higgs scalar remains light, and its properties are precisely 
those of the (weakly-coupled) Standard Model Higgs boson. That is, ~ 0°, with 
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ruho ~ 0{mz), and all other non-minimal Higgs states are significantly heavier than 
rriho. Squared-mass splittings of the heavy Higgs states are of 0{m\)^ which means 
that all heavy Higgs states are approximately degenerate, with mass differences of 
order m\ I "^a^ (here m^o is approximately equal to the common heavy Higgs mass 
scale). In contrast, if the non- minimal Higgs sector is weakly coupled but far from 
the decoupling limit, then is not separated in mass from the other Higgs states. 
In this case, the properties]^ of differ significantly from those of 0°. 

Below, I exhibit the decoupling limit of the most general CP-even two-Higgs- 



doublet model [|13[. It is convenient to define four squared mass combinations: 

ml = Mil cos^ P + -^22 sin^ P + -^12 sin 2/3 , 

1/2 



ml ^ (Ml.Ml, - Mt, 



ml = Ml^ + Ml^ 



m 



g = m^o + m^ , (17) 

in terms of the elements of the neutral CP-even Higgs squared-mass matrix [eq. 
In terms of the above quantities. 



m2 - 1 



ml ± — 4m^o^| — 4m|, 



and 

cos2(/5-a) = ^i^^. (19) 

In the decoupling limit, all the Higgs self-coupling constants Aj are held fixed such 
that Xi ^ 1, while taking m\o ^ Xiv"^. Then ~ 0(f^), and it follows that: 

mfio ~ , mno — m^o ~ mH± , (20) 

and 

111 1 \ 4 

cos'(/5 - a) ^ -J^-L — D _ ^21) 



m 



AO 



Note that eq. (|2T|) implies that cos(/3 — a) = C(m|/m^o) in the decouphng limit, 
which means that the }^ couplings to Standard Model particles match precisely 
those of the Standard Model Higgs boson. These results are easily confirmed by 
considering the cos(/? — a) — > limit of eqs. (P)-([l5|). 

^The basic property of the Higgs coupling strength proportional to mass is maintained. But, 
the precise coupling strength patterns of /i" will differ from those of 0*^ in the non-decoupling 
limit. 
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Although no experimental evidence for the Higgs boson yet exists, there are 
some experimental as well as theoretical constraints on the parameters of the two- 
Higgs doublet model. Experimental limits on the charged and neutral Higgs masses 
have been obtained at LEP. For the charged Higgs boson, m//± > 44 GeV |]14|- 
This is the most model-independent bound and assumes only that the decays 
dominantly into t^i^t, cs and cb. The LEP limits on the masses of and 
are obtained by searching simultaneously for e^e~ —>■ h^ff and e'^e~ h^A^, 
which are mediated by s-channel Z-exchange [15|. The ZZh^ and Zh^A^ couplings 
that govern these two decay rates are proportional to sin(/3 — a) and cos(/5 — a), 
respectively. Thus, one can use the LEP data to deduce limits on rriho and m^o 
as a function of sin(/3 — a). Stronger limits can be obtained in the MSSM where 
sin(/3— a) is determined by other model parameters. At present, taking into account 
data from LEP-1 and the most recent LEP-2 data (at y/s = 161 and 172 GeV), 
one can exclude the MSSM Higgs mass ranges: rriho < 62.5 GeV (independent of 
the value of tan/5) and m^o < 62.5 GeV (assuming tan/3 > 1) |I6[| . 

The experimental information on the parameter tan/? is quite meager. For 
definiteness, let us assume that the Higgs-fermion couplings are specified as in 
eq. (^). The Higgs coupling to top quarks is proportional to gmt/2mw, and is 
therefore the strongest of all Higgs-fermion couplings. For tan/? < 1, the Higgs 
couplings to top-quarks are further enhanced by a factor of 1/ tan/?. As a result, 
some experimental limits on tan/3 exist based on the non-observation of virtual 
effects involving the H~tb coupling. Clearly, such limits depend both on 'mfj± and 
tan/?. The most sensitive limits are obtained from the measurements of B^-B^ 



mixing and the widths of 6 — > 37 and Z ^ bb HT^]. For example, the process 
b s'-f can be significantly enhanced due to charged Higgs boson exchange. If 
there are no other competing non-Standard Model contributions (and this is a big 



if), then present data excludes charged Higgs masses less than about 250 GeV 
(independent of the value of tan/3). In some regions of tan/?, the limits on the 
charged Higgs mass can be even more severe. However, other virtual contributions 
may exist that can cancel the effects of the charged Higgs exchange. For example, in 
the MSSM, constraints on tan/? and m/^± are significantly weaker. For tan/? ^ 1, 
the Higgs couplings to bottom-quarks are enhanced by a factor of tan/?. In this 
case, the measured rate for the inclusive decay oi B ^ X+TPr can be used to set an 
upper limit on tan /3 as a function of the charged Higgs mass. This is accomplished 
by setting a limit on the contribution of the tree-level charged Higgs exchange. 
Present data can be used to set a 2a upper bound of tan/3 < 42(m//±/miy) |TU[ . 
In the MSSM, this bound could be weakened due to one-loop QCD corrections 
mediated by the exchange of supersymmetric particles pO| . 

Theoretical considerations also lead to bounds on tan/3. The crudest bounds 
arise from unitarity constraints. If tan /3 becomes too small, then the Higgs cou- 
pling to top quarks becomes too strong. In this case, the tree-unitarity of processes 
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involving the Higgs-top quark Yukawa coupling is violated. Perhaps this should 
not be regarded as a theoretical defect, although it does render any perturbative 
analysis unreliable. A rough lower bound advocated by Ref. PT|, tan/3 ^ 0.3, 



corresponds to a Higgs-top quark coupling in the perturbative region. A similar 
argument involving the Higgs-bottom quark coupling would yield tan/5 ^ 120. A 
more solid theoretical constraint is based on the requirement that Higgs-fermion 
Yukawa couplings remain finite when running from the electroweak scale to some 
large energy scale A. Above A, one assumes that new physics enters. The limits on 
tan /3 depend on rrit and the choice of the high energy scale A. Using the renormal- 
ization group equations given in Appendix B, one integrates from the electroweak 
scale to A (allowing for the possible existence of a supersymmetry-breaking scale, 
mz < ^susY — determines the region of tan/5-mj parameter space in 

which the Higgs-fermion Yukawa couplings remain finite. This exercise has re- 
cently been carried out at two-loops in Ref. [^. Suppose that the low-energy 
theory at the electroweak scale is the MSSM, and that there is no additional 
new physics below the grand unification scale of A = 2 x 10^^ GeV. Then, for 
rrit = 170 GeV, the Higgs-fermion Yukawa couplings remain finite at all energy 
scales below A if 1.5 ^ tan/3 ^ 65. Note that this result is consistent with the 
scenario of radiative electroweak symmetry breaking in low-energy supersymmetry 
based on supergravity, which requires that 1 ^ tan/? ^ rrit/mb. 



3 The Higgs Sector of the MSSM at Tree Level 

The Higgs sector of the MSSM is a CP-conserving two-Higgs-doublet model, 
with a Higgs potential whose dimension-four terms respect supersymmetry and 
with restricted Higgs-fermion couplings in which $i couples exclusively to down- 
type fermions while $2 couples exclusively to up-type fermions |p. Using the 
notation of eq. (|I]), the quartic couplings Aj are given by 

X, = X, = Xr = 0. (22) 
Inserting these results into eqs. @ and (||), it follows that 

m\o = m^2(tan /5 + cot /3) , 

m|-± = 171%, + , (23) 
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and the tree-level neutral CP-even mass matrix is given by 

^2 _ f rn^o siii^ (3 + cos^ /3 — (m^o + m|) sin/5cos/3 \ 
° y — (m^o + m|) sin P cos P m\o cos^ P + m| sin^ P ) 

The eigenvalues of AIq ^'^^ the squared masses of the two CP-even Higgs scalars 

"^lfo,fto = i ("^Ao + m| ± ^ (m^o + - Am\m\o cos^ 2/? ) . (25) 
and the diagonalizing angle is a, with 



t^o — m| 

2 2 I 1 ^ ^ 



cos2a^-cos2/j( ""f , sin 2a ^ - sin 2/? pf + ) . (26) 



From these results, it is easy to obtain: 



cos\P -a)= . (27) 



Thus, in the MSSM, two parameters (conveniently chosen to be m^o and tan/3) 
suffice to fix all other tree-level Higgs sector parameters. 

Consider the decoupling limit where m^o ^ rriz- Then, the above formulae 
yield 

m\a ~ m| cos^ 2/5 , 

m^o — m\o -\- m| sin^ 2/? , 

2 2,2 

■mjj± = m^o m^y , 

.4 ^;^2 



cos^(/3-a)- , 4 ■ (28) 
4m 



2/Q \ ?^isin 4/3 

Two consequences are immediately apparent. First, m^o ~ m^/o ~ m//±, up to cor- 
rections of C(m|/m^o). Second, cos(/3 — a) = up to corrections of O{m\/m\o). 
Of course, these results were expected based on the discussion of the decoupling 
limit in the general two-Higgs-doublet model given in section 2. 

Finally, a number of important mass inequalities can be derived from the ex- 
pressions for the tree-level Higgs masses obtained above, 

m/jO < m^o 

m/jO < m| cos2/3| < mz , with m = min{mz-i m^o) 
rriHO > mz , 

TTlH^ > TTlw ■ (29) 
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4 Radiative Corrections to the MSSM Higgs Masses 



4.1 Overview 

The tree-level results of the previous section are modified when radiative cor- 
rections are incorporated. Naively, one might expect radiative corrections to have 
a minor effect on the phenomenological implications of the model. However, in the 
MSSM, some of the tree-level Higgs mass relations may be significantly changed 
at one-loop, with profound implications for the phenomenology. For example, 
consider the tree-level bound on the lightest CP-even Higgs boson of the MSSM: 
TTT'hO ^ ''TT'zl cos2/9| < m^. The LEP-2 collider (running at its projected maximum 
center-of-mass energy of 192 GeV, with an integrated luminosity of 150 pb~ ) will 



discover at least one Higgs boson of the MSSM if m^o < niz |15|. Thus, if the tree- 
level Higgs mass bound holds, then the absence of a Higgs discovery at LEP would 
rule out the MSSM. However, when radiative corrections are included, the light 
Higgs mass upper bound may be increased significantly. In the one-loop leading 
logarithmic approximation [ p3| , p^ ] 



m.o ^ rriz cos p + —5 — 5- In — , (30) 

where M^^, Mi^ are the masses of the two top-squark mass eigenstates. Observe 
that the Higgs mass upper bound is very sensitive to the top mass and depends 
logarithmically on the top-squark masses. In addition, due to the increased upper 
bound for m/jo , the non-observation of a Higgs boson at LEP-2 cannot rule out the 
MSSM. 

Although eq. (|30| ) provides a rough guide to the Higgs mass upper bound, 
it is not sufficiently precise for LEP-2 phenomenology, whose Higgs mass reach 
depends delicately on the MSSM parameters. In addition, in order to perform 
precision Higgs measurements and make comparisons with theory, more accurate 
results for the Higgs sector masses (and couplings) are required. The radiative 
corrections to the Higgs mass have been computed by a number of techniques, and 
using a variety of approximations such as effective potential [ pi| , p5| , pBtP7| , P^ ] and 



diagrammatic methods |23,29,30,31,32,33,34|. Complete one- loop diagrammatic 



computations of the MSSM Higgs masses have been presented by a number of 



groups the resulting expressions are quite complex, and depend on all the 
parameters of the MSSM. (The dominant two-loop next-to-leading logarithmic 
results are also known [Q.) Moreover, as noted above, the largest contribution 
to the one-loop radiative corrections is enhanced by a factor of ml and grows 
logarithmically with the top squark mass. Thus, higher order radiative corrections 
can be non-negligible for large top squark masses, in which case the large logarithms 
must be resummed. 
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The renormalization group (RG) techniques for resumming the leading loga- 



rithms has been developed by a number of authors ||35| . |36| . |37[| . The computation 
of the RG-improved one-loop corrections requires numerical integration of a cou- 
pled set of RG equations |^ . Although this program has been carried out in the 



literature, the procedure is unwieldy and not easily amenable to large-scale Monte- 
Carlo analyses. Recently, two groups have presented a simple analytic procedure for 
accurately approximating m/jO. These methods can be easily implemented, and in- 
corporate both the leading one-loop and two-loop effects and the RG-improvement. 
Also included are the leading effects at one loop of the supersymmetric thresholds 
(the most important effects of this type are squark mixing effects in the third gen- 
eration). Details of the techniques can be found in Refs. |^ and p8|. Here, I 



simply quote two specific bounds, assuming rrit = 175 GeV and ^ 1 TeV: 
m/jO ^ 112 GeV if top-squark mixing is negligible, while m/jO ^ 125 GeV if 
top-squark mixing is "maximal" . Maximal mixing corresponds to an off-diagonal 
squark squared- mass that produces the largest value of niho. This mixing leads to 
an extremely large splitting of top-squark mass eigenstates. 

The charged Higgs mass is also constrained in the MSSM. At tree level, mH± 
is given by eq. (^31), which implies that charged Higgs bosons cannot be pair pro- 
duced at LEP-2. Radiative corrections modify the tree- level prediction, but the 
corrections are typically smaller than the neutral Higgs mass corrections discussed 
above. Although m/^± > mw is not a strict bound when one-loop corrections are 
included, the bound holds approximately over most of MSSM parameter space (and 
can be significantly violated only when tan/5 is well below 1, a region of parameter 
space that is theoretically disfavored). 

In the remainder of this section, I shall present formulae which exhibit the 
leading contributions to the one-loop corrected Higgs masses. Symbolically, 

M' = Ml^^ + AMl^, (31) 

where the subscript ILL refers to the tree-level plus the one-loop leading logarith- 
mic approximation to the full one-loop calculation, and the subscript mix refers to 
the contributions arising from qi-qR mixing effects of the third generation squarks. 
The CP-even Higgs mass-squared eigenvalues are then obtained by using eq. (|^ 
and the corresponding mixing angle, a, is obtained from eq. (§). 

In the simplest approximation, squark mixing effects are neglected and the su- 
persymmetric spectrum is characterized by one scale, called Mg^gy- We assume 
that MgusY is sufficiently large compared to mz such that logarithmically enhanced 
terms at one-loop dominate over the non-logarithmic terms.Q In this case, the 

^If this condition does not hold, then the radiative corrections would constitute only a minor 
perturbation on the tree-level predictions. 
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full one-loop corrections {e.g., obtained by a diagrammatic computation) are well 
approximated by the one-loop leading logarithmic approximation. Next, we incor- 
porate the effects of squark mixing, which constitute the largest potential source 
of non-logarithmic one-loop corrections. In particular, these contributions to the 
Higgs mass radiative corrections arise from the exchange of the third generation 
squarks. Now, the approximation is parameterized by four supersymmetric param- 
eters: MgygY (a common supersymmetric particle mass) and the third generation 
squark mixing parameters: At, A), and fi. A more comprehensive set of formu- 
lae can be derived by treating the third generation squark sector more precisely 
by accounting for non-degenerate top and bottom squark masses. This approx- 
imation is characterized by seven supersymmetric parameters — the three squark 
mixing parameters mentioned above, three soft-supersymmetry-breaking diagonal 
squark mass parameters, Mq, Mu, and Mjj, and a common supersymmetry mass 
parameter Mg^gY which characterizes the masses of the first two generations of 
squarks, the sleptons, the charginos, and the neutralinos. 

Given an approximation to the one-loop Higgs mass (at some level of ap- 
proximation as described above), one must incorporate the RG-improvement if 
MgugY 3> mz. A simple analytic procedure of Ref. is described in the sec- 
tion 5, and some numerical results are presented there. Similar results have also 
been obtained by Carena and collaborators, where analytic approximations to the 
RG-improved radiatively corrected MSSM Higgs masses are also developed p8 



Although the approaches are somewhat different, the numerical results (in cases 
which have been compared) typically agree to within 1 GeV in the evaluation of 
Higgs masses. 



4.2 One-Loop Leading Logarithmic Corrections to the MSSM Higgs 
Masses 

The leading logarithmic expressions for Higgs masses can be computed from 
the one-loop renormalization group equations (RGEs) of the gauge and Higgs self- 
couplings, following Ref. [^. The method employs eqs. (^) and (^), which are 
evaluated by treating the Aj as running parameters evaluated at the electroweak 
scale, M^gg^]^. In addition, we identify the W and Z masses by 



2 12/2, 2\ 

= a9 (^1 +^^2) . 



(32) 



where the running gauge couplings are also evaluated at M^^^^. Of course, the 
gauge couplings, g and g' are known from experimental measurements which are 
performed at the scale M^^^^. The Aj(M^^^j^) are determined from supersymmetric 
boundary conditions at Mg^gY and RGE running down to M^^^-^. That is, if 
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super symmetry were unbroken, then the Aj would be fixed according to eq. (0). 
Since supersymmetry is broken, we regard eq. (^) as boundary conditions for the 
running parameters, vahd at (and above) the energy scale Afg^gY- That is, we take 

^i(^susy) = ^^(MI^sy) = i[^'(^susY) +^''(^susy)]' 
A3(Mg^,,g^) = i f/(M|^gY) - ^?'^(Mg^^g^, 



^4(M|uSy) = -i^^'(^SUSY)' 

^slMg^usY) = >^(^(Msvsy) = ^^iM^vsy) = ' (33) 

in accordance with the tree-level relations of the MSSM. At scales below Mg^gy, 
the gauge and quartic couplings evolve according to the renormalization group 
equations (RGEs) of the non-supersymmetric two-Higgs-doublet model given in 
eqs. (B.5)-(B.7). These equations are of the form: 

-j^ = Pi{pi,P2, . . .) witht = ln/i^, (34) 

where fi is the energy scale, and the Pi are the parameters of the theory {pi = 
g^, \k, ■ ■ ■)■ The relevant /3-functions can be found in Appendix B. The boundary 
conditions together with the RGEs imply that, at the leading-log level, A5, Xq 
and Ay are zero at all energy scales. Solving the RGEs with the supersymmetric 
boundary conditions at Mg^gy, one can determine the Aj at the weak scale. The 
resulting values for K^M^^^^J are then inserted into eqs. @ and (H) to obtain the 
radiatively corrected Higgs masses. Having solved the one-loop RGEs, the Higgs 
masses thus obtained include the leading logarithmic radiative corrections summed 
to all orders in perturbation theory. 

The RGEs can be solved by numerical analysis on the computer. In order 
to derive the one-loop leading logarithmic corrections, it is sufficient to solve the 
RGEs iteratively. In first approximation, we can take the right hand side of eq. ( p^ 
to be independent of /i^. That is, we compute the Pi by evaluating the parameters 
Pi at the scale n = Mg^gy. Then, integration of the RGEs is trivial, and we obtain 

MKeJ = P^iM^vsr) - A In (^) . (35) 

\ weak / 

This result demonstrates that the first iteration corresponds to computing the one- 
loop radiative corrections in which only terms proportional to InM^^^^ are kept. 
It is straightforward to work out the one-loop leading logarithmic expressions for 
the Aj and the Higgs masses. 
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First consider the charged Higgs mass. Since A5(/i^) = at all scales, we need 
only consider A4. Evaluating at /x = Mg^gY, "we compute 



A4(m^) 



327r2 



4jV 



9' + 59^9" 



^9' 



+ 



22 4 



(36) 



w 



The terms proportional to the number of generations A^^ = 3 and the number of 
Higgs doublets Nh = 2 that remain in the low-energy effective theory at the scale 
H = m\Y have their origin in the running of g'^ from Mg^gY down to m\Y. In deriving 
this expression, I have taken M^^^^ = mw This is a somewhat arbitrary decision, 
since another reasonable choice would yield a result that differs from eq. (|36D by 
a non-leading logarithmic term. Comparisons with a more complete calculation 
show that one should choose M^^^^. = rriw in computations involving the charged 
Higgs (and gauge) sector, and M^^^^^ = rriz in computations involving the neutral 
sector. 

The above analysis also assumes that rrit ~ 0{my/). Although this is a good 
assumption, we can improve the above result somewhat by decoupling the (t, h) 
weak doublet from the low-energy theory for scales below mt . The terms in eq. ( ^61) 
that are proportional to and/or ml arise from self-energy diagrams containing 
a th loop. Thus, such a term should not be present for my/ < n <mt. In addition, 
we recognize the term in eq. ( p6| ) proportional to the number of generations Ng as 
arising from the contributions to the self-energy diagrams containing either quark 
or lepton loops (and their supersymmetric partners). To identify the contribution 
of the th loop to this term, simply write 



Ng = \Ng{N, + 1) = \N, + \[N,{Ng " 1 ) + A^, 



a\ 1 



(37) 



where A^c = 3 colors. Thus, we identify ^N^. as the piece of the term proportional 
to Ng that is due to the th loop. The rest of this term is then attributed to the 
lighter quarks and leptons. Finally, the remaining terms in eq. (|3^) are due to the 



contributions from the gauge and Higgs boson sector. The final result is [31 



A4(m^) 



12 N,9' 
2" 



327r2 



mt 



2m^ 



+ 



2 2 
m^m^ 



s}c}m' 



mt 



967r2 



{ [N,{Ng -1) + Ng+ '^Nh - 10] + 15(7 V'} In 



'M2 ^ 

SUSY 



m 



w 



(38) 
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Inserting this result (and A5 = 0) into eq. (H), we obtain the one- loop leading- 
logarithmic (ILL) formula for the charged Higgs mass 



Ncg' 



327r2m^ 



2 ''^W 




X In 



SUSY 



+ 



w 



{ [iV,(iV, _ 1) + iv, - 9] (7' + 15(?''} In 



SUSY 



(39) 



Since this derivation makes use of the two-Higgs-doublet RGEs for the Aj, there 
is an implicit assumption that the full two-doublet Higgs spectrum survives in the 
low-energy effective theory at yU = my/. Thus, I have set = 2 in obtaining 
eq. (^) above. It also means that m^o cannot be much larger than mi4/.0 The 
leading logarithms of eq. (PPf) can be resummed to all orders of perturbation theory 



by using the full RGE solution to \^{m^^ 



(40) 



Although the one-loop leading-log formula for mii± [eq. (|39|)] gives a useful 
indication as to the size of the radiative corrections, non-leading logarithmic con- 
tributions can also be important in certain regions of parameter space. A more 



complete set of radiative corrections can be found in the literature [26,29,30,31,34 



However, it should be emphasized that the radiative corrections to the charged 
Higgs mass are significant only for tan/3 < 1, a region of MSSM parameter space 
not favored in supersymmetric models. 

The computation of the neutral CP-even Higgs masses follows a similar proce- 
dure. The results of Ref. |3^ are summarized below. From eq. (|^), we see that we 
only need results for Ai, A2 and A3 = A3 + A4 + A5. (Recall that A5 = Ag = A7 = 
at all energy scales.) By iterating the corresponding RGEs as before, we find 

and the 
For 



O(MgugY), then and would aU have masses of order A/gygy, 

effective low-energy theory below -Mgygy would be that of the minimal Standard Model, 
example, for m^o = -^susy' leading logarithmic corrections to the charged Higgs mass can 



be obtained from to^^^ = m^o + by treating as a running parameter evaluated at m^o . 



Re-expressing tovi/(toao) in terms of the physical W mass yields the correct one- loop leading 
log correction to m'^^±. For mz < to^o < Afgygy, one can interpolate between the effective 
two-Higgs doublet model and the effective one-Higgs doublet model. 
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Ai(m|) = 



A2(m|) 



A3(m|) 




+ f 12iV,^J^ - 6iV,^J^ + + p + p + P2^ 




+ fl2iV,^-6^„^ + P,)luf^)' 



In 



(41) 



where 



Pt 
Pf 



2H 



P' 



K{N, - 1)[2 - + S{el + e^)s4^] + 7V,[2 - + 8s^^] , 



-44 + 106s^ - 62s^ , 
10 + 34s^-26s^, 
-10 + 2s^-2s^, 



= 8 - 22s^ + 10s 



(42) 



In the above formulae, the electric charges of the quarks are Ct = 2/3, = —1/3, 
and the subscripts t, b, f, g and 2H indicate that these are the contributions from 
the top and bottom quarks, the other fermions (leptons and the first two genera- 
tions of quarks), the gauge bosons and the Higgs doublets, and the corresponding 
supersymmetric partners, respectively. 

As in the derivation of A4(m^) above, we have improved our analysis by remov- 
ing the effects of top-quark loops below fi = rrit. This requires a careful treatment 
of the evolution of g and g' at scales below ji — rrit. The correct procedure is 
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somewhat subtle, since the full electroweak gauge symmetry is broken below top- 
quark threshold; for further details, see Ref . . However, the following pedestrian 
technique works: consider the RGE for + g'"^ valid for n < Mg^gy 



d 
di 



+ 9'' 



1 



(43) 



This equation is used to run g'^ + g'^, which appears in eq. (^3]), from Mg^gy down 
to rriz- As before, we identify the term proportional to A^^^ as corresponding to the 
fermion loops. We can explicitly extract the t-quark contribution by noting that 



N, {8g' + fg" 



9'N, 



64^4 _1Ro2 I 
3 ^'^^W ' 



N,[l + {Ng - 1)](1 - 4eiS^ + %e\s%, 



+ N,N,{\ + 4ebS^ + %eis%,) + Ar,(2 - 4s^ + 8s 



W) ( ' 



(44) 



where in the first line of the last expression, the term proportional to 1 corresponds 
to the t-quark contribution while the term proportional to A^^^ — 1 accounts for the 
u and c-quarks; the second line contains the contributions from the down-type 
quarks and leptons respectively. Thus, iterating to one-loop. 



[g' + g'^mi^^^) = [g' + g'']im',) + 



967r2c^ 



PAn 



+ 



A + Py + (4 + 4)Ar^-44c: 



In 



SUSY 



(45) 



Again, we take A^^ = 2, since the low-energy effective theory between mz and 
MgusY consists of the full two-Higgs doublet model. Eq. (^5]) was used in the 
derivation of eq. (^Tf). 

We now return to the computation of the one-loop leading log neutral CP-even 
Higgs squared-mass matrix. The final step is to insert the expressions obtained in 
eq. (|4l|) into eq. (^. The resulting matrix elements for the mass-squared matrix 
to one-loop leading logarithmic accuracy are given by 



18 



f Ka2 \ 2 2, 2 2, d^^'z^^P 



In 



V ml ) 




+ f 12iV,^J^ - 6iV,^^ + + p + p + p^^ ) In 



/ XA'i \ 2 2, 2 2, 




(A<i2)iLL = -spcpi^ m^o + + ^^^2(0^ 



Pt-3iV,^ 




+ ( - 3iV,^J^ + p, + p + p; + p'^ In 




(46) 



Diagonalizing this matrix [eq. (^6[) ] yields the radiatively corrected CP-even Higgs 
masses and mixing angle a. 

The analysis presented above assumes that m^o is not much larger than 0{mz) 
so that the Higgs sector of the low-energy effective theory contains the full two- 
Higgs-doublet spectrum. On the other hand, if m^o ^ rriz, then only remains in 
the low-energy theory. In this case, we must integrate out the heavy Higgs doublet, 
in which case one of the mass eigenvalues of A^q bl- (2D] is much larger than the 
weak scale. In order to obtain the effective Lagrangian at M^^^^^, we first have to 
run the various coupling constants to the threshold rriAO- Then we diagonalize the 
Higgs mass matrix and express the Lagrangian in terms of the mass eigenstates. 
Notice that in this case the mass eigenstate is directly related to the field with the 
non-zero vacuum expectation value [i.e., /3(m^o) = a{mAo) + tt/2 + O{m%/m\o)]. 

Below m^o only the Standard Model Higgs doublet (p = C/3$i + Sj3^2 remains. 
The scalar potential is 



V = mJ(0V) + iA(0V)2, 

and the light CP-even Higgs mass is obtained using m^o 
Standard Model for A is 



(47) 

The RGE in the 



6A^ + | 



2g* + {g^ + g 



i2\2 



X 



1/ + 



2YNr.h^i 



(48) 
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where the summation is over all fermions with hj^ = gm fj (\/2mw) ■ The RGEs 
for the gauge couplings are obtained from j3g2 and /3g'2 given in Appendix B by 
putting Nh = 1. In addition, we require the boundary condition for A at m^o 

A(mAo) = c^Ai + s^A2 + 2s^c^(A3 + A4 + As) + 4c^s^A6 + 4c;3S^A7 {niAo) 



9^ 



3847r^cV 



In 



SUSY 



m 



AO 



X 



ml 



mr 



mt 



mt 



rrir. 



nir 



+ clp{P, + A + Pf) + (4 + 4)(P, + P,h) - '^slc%P', + P2h) 



, (49) 



where {g + 9 )clf^ is to be evaluated at the scale m^o as indicated. The RGE for 
g^ + g'^ was given in eq. (^); note that at scales below we must set Nh = 1- 
Finally, we must deal with implicit scale dependence of C2^. Since the fields $j 
[i = 1, 2) change with the scale, it follows that tan (3 scales like the ratio of the two 
Higgs doublet fields, i.e., 

citan^/5 ^jd 



(50) 



tan^ f5 dt 

Thus we arrive at the RGE for cos 2/? in terms of the anomalous dimensions 7, 



given in eq. (p.6| ). Solving this equation iteratively to first order yields 



(niAo) = C2p{mz) + 4c2/3C^s^(7i - 72) In 



AO 



mr 



(51) 



The one loop leading log expression for m\o = \{mz)v'^ can now be obtained 
by solving the RGEs above for \{mz) iteratively to first order using the boundary 
condition given in eq. (|49|). The result is 



("^feo)lLL 



g^ml 



+ 



ml 



P, + P2h) (4 + 4) - 2.314 [P9 + P2H ) In 
ml „ , ^ m 



6N,C2f3^ + 4f,{Pb + Pf) 



4 2 
12iVe^ + 6N,C2(3^ + clpPt 



mr 



mr 



In 



M2 
SUSY 

2 

mf 



{4 + 4) 



P-2H - 2CijSpP2H 



Pi 



IH 



~ 

SUSY 

ml 




(52) 
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where the term proportional to 



P,H = -94^ + (1 - 24 + 2sty)c% , (53) 

corresponds to the Higgs boson contribution in the onc-Higgs-doublet modeL The 
term in eq. (0) proportional to In(m^o) accounts for the fact that there are two 
Higgs doublets present at a scale above m^o but only one Higgs doublet below 

777.^0. 

We can improve the above one-loop leading log formulae by reinterpreting the 
meaning of Mg^gy. For example, all terms proportional to ln{M^^^-y / mf) arise 
from diagrams with loops involving the top quark and top-squarks. Explicit di- 
agrammatic computations then show that we can reinterpret M^^^^ = Mi^Mi^. 
Note that with this reinterpretation of Mg^gy, the top quark and top squark loop 
contributions to the Higgs masses cancel exactly when Mi^ = Mf^ = rrit, as required 
in the supersymmetric limit. Likewise, in terms proportional to or powers of 
rUb multiplied by \iii{M^^^y / m"^^) , we may reinterpret Mg^gy = M^^Mi^. Terms 
proportional to P f In^M"^^^^ / im?z) come from loops of lighter quarks and leptons 
(and their supersymmetric partners) in an obvious way, and the corresponding 
Mgygy can be reinterpreted accordingly. The remaining leading logarithmic terms 
arise from gauge and Higgs boson loops and their supersymmetric partners. The 
best we can do in the above formulae is to interpret Mg^gy as an average neu- 
tralino and chargino mass. To incorporate thresholds more precisely requires a 
more complicated version of eq. (0), which can be easily derived from formulae 
given in Ref. [E^. The explicit form of these threshold corrections can be found 



in Ref. . However, the impact of these corrections are no more important than 
the non-leading logarithmic terms which have been discarded. 

The largest of the non- leading logarithmic terms is of 0{g'^mf), which can 
be identified from a full one-loop computation as being the subdominant term 
relative to the leading 0{g'^m^ In M^^^y) term in A^22- Thus, we can make a 
minor improvement on our computation of the one-loop leading-log CP-even Higgs 
squared mass matrix by taking 



where A^^ll is the matrix whose elements are given in eq. (^). One can check 
that this yields at most a 1 GeV shift in the computed Higgs masses. 

4.3 Leading Squark Mixing Corrections to the MSSM Higgs Masses 

In the case of multiple and widely separated supersymmetric particle thresh- 
olds and/or large squark mixing (which is most likely in the top squark sector). 
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new non-leading logarithmic contributions to the scalar mass-squared matrix can 



become important. As shown in Ref. ^7j, such effects can be talcen into account 



by modifying the boundary conditions of the A, at the supersymmetry breaking 
scale [eq. ([33|)], and by modifying the RGEs to account for multiple thresholds. In 
particular, we find that A5, Xq and A7 are no longer zero. If the new RGEs are 
solved iteratively to one loop, then the effects of the new boundary conditions are 
simply additive. 

In this section, we focus on the effects arising from the mass splittings and g^- 
mixing in the third generation squark sector. The latter generates additional 
squared-mass shifts proportional to and thus can have a significant impact on 
the radiatively corrected Higgs masses |]2^. First, we define our notation (we follow 
the conventions of Ref. 0). In third family notation, the squark mass eigenstates 
are obtained by diagonalizing the following two 2x2 matrices. The top-squark 
squared-masses are eigenvalues of 



Mq + ml + tLVfil mtXt 



(55) 



\^ mtXt Ml + ml + tRml J ' 

where Xt = At — /icot/3, = — ej sin^ 6*14/) cos 2/5 and tR = sin^ 6^1^ cos 2/?. 
The bottom-squark squared-masses are eigenvalues of 

/ Mi^ + ml + bLml ^bXh \ 

[ mbXb Ml + ml + hRml ) ' 



(56) 



where Xf, = Af, — ^ tan /5, 6^ = (— | — e?, sin^ 9w) cos 2/? and Br = Cb sin^ 9w cos 2/?. 
Mq, Mu, Md, At, and Ab are soft-supersymmetry-breaking parameters, and /i is 
the supersymmetric Higgs mass parameter. We treat the squark mixing pertur- 
batively, assuming that the off-diagonal mixing terms are small compared to the 
diagonal terms. 

At one-loop, the effect of the squark mixing is to introduce the shifts AA^^j^ 
and (Am^±)j^j^. In order to keep the formulae simple, we take Mq = Mu = 
Md = ^susY' where Mg^gY is assumed to be large compared to niz- Thus, the 
radiatively corrected Higgs mass is determined by m^o, tan/3, Mg^gy, At, Ab, and 
/X. The more complex case of non- universal squark squared-masses (in which Mq, 



Mu, and Md are unequal but still large compared to mz) is treated in Ref. ^8 
It is convenient to define 

Xt = At — fi cot P , Yt = At + fi tan /? , 

Xb = Ab — ^tan P , Yb = Ab + fi cot P . (57) 

We assume that the mixing terms rritXt and rribXb are not too large.0 Then, the 

^Formally, the expressions given in eqs. (^8|)-(|6^) are the results of an expansion in the variable 
[Ml — M|)/(Af|^ + Af|), where Mf, M| are the squared-mass eigenvalues of the squark mass 
matrix. Thus, we demand that mtXt/Mg^^-y ^ 1. For example, for Mgygy = 1 TeV, values of 
Xt/MgygY ^ 3 should yield an acceptable approximation based on the formulae presented here. 
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elements of the CP-even Higgs squared-mass matrix are given by: 

= Ml^^ + AMI 
where A^ill has been given in eq. (^6]), and 



(2. 

mix ' 



(5^ 



{AMI, 



1 



12M2 



cos^ /? 

m^mlAiy^Xi, + \Aij) — m\mljji cot l3{Xt + |/i cot /3) 



^mfAtXt 
sin^ P 



AtXt 



12M2 



-m\mlAt{Xt + i^t) — m|m^/itan/?(X{, + |/itan/?) 



3M|ugY sin2/3 



3MI^^^ cos2 



(AM 



4m^/xX, / _ AXt \ AmtfiX, / _ A,X, 



sin^ /? \^ 6M|ygY y ~'~ cos^ /? 



^^^^SUSY. 



-m|m^ cot /5 XtYt + ^(/i^ + A^^) - m|m^ tan P XhYf, + ^(/i^ + A^' 



• (59) 



If mz -C m^o < MgugY) ^ separate analysis is required. One finds that eq. 
is shifted by 



(Am 



2m^X2 ( 1 ^ 



+ 2m^xMl 



Xi 



SUSY- 



12M2 



+ |m| cos 213 [m^t (x^^ + 1(A^ -/x^ cot^ - (x^^ + |(A^-/i^ tan^ I . (60) 



Squark mixing effects also lead to modifications of the charged Higgs squared- 
mass. One finds that the charged Higgs squared-mass obtained in eq. ( pQ]) is shifted 
by 



2 ^ N,g' 



2mfmlr{fi^-2Al) , 2mlmlr{fi^ - 2A 



sin^ p 



cos^ /? 



— 3/i 
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mt 
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SUSY 



(61) 
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5 RG-Improvement and Numerical Results for the MSSM Higgs Masses 



The RG-improved Higgs masses (in the absence of squark mixing) are com- 
puted by solving the set of coupled REGs for the Aj(M^^^j^), subject to the bound- 
ary conditions specified in eq. (|3^). Squark mixing effects are incorporated into 
the procedure by modifying the boundary conditions as described in Ref. ^7\. 
Hempfling, Hoang and 1 found a simple analytic algorithm which reproduces 
quite accurately the results of the numerical integration of the RGEs. 

The procedure starts with the formulae of section 4. The Higgs masses take 
the form given symbolically in eq. (pi]). Then, 



M 



IRG 



ILL 



-^ILL [mt{fXt),mb{fXb)] + AA^mix K(/^ 



(62) 



where 



SUSY ' 



SUSY ' 



N = ^susY = i,b). (63) 



That is, the numerically integrated RG-improved CP-even Higgs squared-mass 
matrix, A^irq? ^^^^ approximated by replacing all occurrences of and rrif, in 
M iLhi'^^t, rrib) and AA^^j^(mt, m^) by the corresponding running masses evaluated 
at the scales as indicated above.Q To implement the above algorithm, we need 
formulae for mb{fi) and mt{f^). First, consider uiao = 0{mz)- In this case, at 
mass scales below M, 



SUSY' 



the effective theory of the Higgs sector is that of a 
non-supersymmetric two-Higgs-doublet model. In this model, the quark mass is 
the product of the Higgs-quark Yukawa coupling (hq) and the appropriate Higgs 
vacuum expectation value: 



rribilj) = ■^hb{fi)vi{fi) 



(64) 



At scales /x < Mg^gY, we employ the one- loop non-supersymmetric RGEs of the 
two-Higgs doublet model for hb, hf, and the vacuum expectation values Vi and V2 
(see Appendix B). This yields 
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din /I 



64 7r2 
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Qhl + 2hl - 32(7, 
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(65) 



^In this section, an overline above a quantity will indicate that the replacement of mt and nib 
by the appropriate running mass has been made. 
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For rriAo = C'(^susy)' effective theory of tlie Higgs sector at mass scales 
below MgygY is that of the one-Higgs doublet Standard Model. In this case, we 
define mq{fx) = /i^^(/i)w(/i)/\/2, where v{mz) — 246 GeV is the one-Higgs-doublet 
Standard Model vacuum expectation value. In this case eq. (|65D is modified by 
replacing Qhj + 2hl with Q{hf^Y — Qijif^)'^ in the RGE for (and interchange b 
and t to obtain the RGE for ml). 

To solve these equations, we also need the evolution equations of Qs, and g' . 
But, an approximate solution is sufficient for our purposes. Since g' is small, we 
drop it. We do not neglect the hh dependence which may be significant if tan/3 
is large. Then, we can iteratively solve eq. ( |65D to one loop by ignoring the 
dependence of the right hand side. We find 



1-i 



as - j^{ab + 3at) \n {jj,^ / mf) , m^o ~ 0(mz) 



(66) 

where at = hf/4:7!-, etc., and all coupling on the right hand side are evaluated at 
rrit. Similarly, 



1-i 



^{at + 3ab) \ ln(/iV"i|) > i^ao - 0{mz), 
ln(/iVm|) , ~ C(Mg 



iQ {af af 



'SUSY/'' 

(67) 

For intermediate values of m^o, one may interpolate the above formulae between 
the two regions. Using eqs. (|66|) and (|67D in eq. (|62D , and diagonalizing the result- 
ing squared-mass matrix yields our approximation to the RG-improved one-loop 
neutral CP-even Higgs squared- masses. 

We may also apply our algorithm to the radiatively corrected charged Higgs 
mass. However, in contrast to the one-loop radiatively corrected neutral Higgs 
mass, there are no one-loop leading logarithmic corrections to that are pro- 
portional to mf. Thus, we expect that our charged Higgs mass approximation will 
not be quite as reliable as our neutral Higgs mass approximation. 

Let us now compare various computations of the one-loop corrected light CP- 
even Higgs mass. In the first set of examples, all squark mixing effects are ignored. 
First, we evaluate two expressions for the RG-unimproved one- loop Higgs mass — 
the one- loop leading log Higgs mass calculated from A^^ll ^"^^ from a simplified 
version of which only the dominant terms proportional to are kept. 

In the latter case, we denote the neutral CP-even Higgs squared-mass matrix by 
MluY = Ml + AMIjj^, where 



25 



In many analyses of A^ilt ^"^^ -^ill that have appeared previously in the 
literature, the Higgs mass radiative corrections were evaluated with the pole mass, 
rrit. Some have argued that one should take rrit to be the running mass evaluated 
at rrit, although to one-loop accuracy, the two choices cannot be distinguished. 
Nevertheless, because the leading radiative effect is proportional to mf, the choice 
of rut in the one- loop formulae is numerically significant, and can lead to differences 
as large as 10 GeV in the computed Higgs mass. In Ref. [^, the choice of using 
mt{mt) as opposed to m\°^'^ (prior to RG- improvement) is justified by invoking 
information from a two-loop analysis. Thus, our numerical results for the light 
CP-even Higgs mass before RG-improvement are significantly lower (when Mg^gy 
is large) as compared to the original computations given in the literature, for fixed 
mf°'^. We have taken mt{mt) = 166.5 GeV in all the numerical results exhibited 
below. We then apply our algorithm for RG-improvement by replacing rut and 
by the appropriate running masses as specified in eqs. (p^-(|63D. 

We now show examples for m^o = 1 TeV and two choices of tan/? in Fig. |l| 
[tan/3 = 20] and Fig. | [tan/3 = 1.5], and for m^o = 100 GeV and tan/3 = 20 in 
Fig. ^.[] Each plot displays five predictions for m/jO based on the following methods 
for computing the Higgs squared-mass matrix: (i) A^ilt! (h) AiiLL) (hi) ■M.'^ilt] 
(iv) AI^ill; and (v) Mf^Q [the overline notation is defined in the footnote below 
eq. (0)]. The following general features are noteworthy. First, we observe that over 
the region of Mg^gy shown, A^irg — ■M'^ill- In fact, m/^o computed from A^^ill 
is within 1 GeV of the numerical RG-improved m/jO in all sensible regions of the 
parameter space (1 < tan/3 < rrit /mi, and rrit, < ^susy — 2 TeV). For values 
of Mgygy > 2 TeV, the Higgs masses obtained from AI^ill begin to deviate from 
the numerically integrated RG-improved result. Second, the difference between 
rriho computed from AI^ll ^^d from Ali^c is non- negligible for large values of 
Mgugy; neglecting RG-improvement can lead to an overestimate of m/jO which in 
some areas of parameter space can be as much as 10 GeV. Finally, note that while 
the simplest approximation of rriho based on AIilt reflects the dominant radiative 
corrections, it yields the largest overestimate of the light Higgs boson mass. 

We next consider some examples in which squark-mixing effects are included. 
As above, we compare the value of rrihO computed by different procedures. Prior to 
RG-improvement, we first compute m/jO by diagonalizing AI^ll + ^-^mix- Next, 
we perform RG-improvement as in Ref. by numerically integrating the RGBs 
for the Higgs self-couplings and inserting the results into eq. (^; the resulting CP- 
even scalar squared-mass matrix is denoted by MI^q. Finally, we extract rrihO and 
compare it to the corresponding result obtained by diagonalizing Ai^iLLH- AAi^mix 
given by eq. (|62D . These comparisons are exhibited in a series of figures. First, we 
plot rriho vs. Xi/Mgygy for Mg^gy = rriAo = — /x = 1 TeV for two choices of tan/3 

'''For rrij^o = 100 GeV and tan/3 = 1.5, the resulting light Higgs mass lies below experimental 
Higgs mass bounds obtained by the LEP collaborations ||l^ . 
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Fig. 1. The radiatively corrected light CP-even Higgs mass is plotted 
as a function of Mg^gy for tan /? = 20 and m^o = 1 TeV. The one- loop 
leading logarithmic computation [dashed line] is compared with the RG- 
improved result which was obtained by numerical analysis [solid line] 
and by using the simple analytic result given in eq. (^) [dot-dashed 
line]. For comparison, the results obtained using the leading mf ap- 
proximation of eq. (|6^ ) [higher dotted line], and its RG-improvement 
[lower dotted line] are also exhibited. -^gusY characterizes the scale 
of supersymmetry breaking and can be regarded (approximately) as a 
common supersymmetric scalar mass; squark mixing effects are set to 
zero. The running top quark mass used in our numerical computations 
is mt{mt) = 166.5 GeV. All figures are taken from Ref. 



in Fig. I [tan/5 = 20] and Fig. || [tan/5 = 1.5]. Note that Fig. | is of particular 
interest, since it allows one to read off the maximal values of rriho as a function of 
Xt for MgugY < 1 TeV, which were quoted in section 4.1. The maximum value of 
the Higgs mass occurs for ~ 2.4MgugY. 

The reader may worry that this value is too large in light of our perturbative 
treatment of the squark mixing. However, comparisons with exact diagrammatic 
computations confirm that these results are trustworthy at least up to the point 
where the curves reach their maxima. From a more practical point of view, such 
large values of the mixing are not very natural; they cause tremendous splitting 
in the top-squark mass eigenstates and are close to the region of parameter space 
where the SU(2) xU(l) breaking minimum of the scalar potential becomes unstable 
relative to color and/or electromagnetic breaking vacua [^ . 
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Fig. 2. The radiatively corrected light CP-even Higgs mass is plotted as 
a function of Mg^gy for tan (5 = 1.5 and m^o = 1 TeV. See the caption 
to Fig. |1[ 




Fig. 3. The radiatively corrected light CP-even Higgs mass is plotted 
as a function of Mg^gy for tan /? = 20 and rrij^o = 100 GeV. See the 
caption to Fig. |l[ 
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Fig. 4. The radiatively corrected light CP-even Higgs mass is plotted as 
a function of Xt/Mg^gy, where Xt = At — [i cot /3, for Afg^gy = m^o = 
—[I = 1 TeV and tan/? = 20. See the caption to Fig. ||. 
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Fig. 5. The radiatively corrected light CP-even Higgs mass is plotted as 
a function of Xt/Mg^gyi where Xt = At — fi cot /3, for Mg^gY = "^ao = 
—fj, = 1 TeV and tan/? = 1.5. See the caption to Fig. 0. 
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Fig. 6. The radiatively corrected, RG-improved light CP-even Higgs mass is plotted as 
a function of Xj/Mg^gy, where Xt = At — /icot/3, for Mg^gy = m^o = 1 TeV and 
two choices of tan/3 = 1.5 and 20. Three values of fi are plotted in each case: —1 TeV 
[dashed], [solid] and 1 TeV [dotted]. Here, we have assumed that the diagonal squark 
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Fig. 7. The radiatively corrected, RG-improved light CP-even Higgs mass is plotted as 
a function of X^/Mgusy for Afgugy = 1 TeV and m^o = 100 GeV. See the caption to 
Fig. |. 
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Fig. 8. The radiatively corrected, RG-improved light CP-even Higgs mass is plotted as a 
function of -Mg^gy for Xt = 2.4MgugY for three choices of (tan/3, m^o)= (20,1), (1.5,1), 
and (1.5,0.1), where m^o is specified in TeV units. The solid line depicts the numerically 



integrated result, and the dot-dashed line indicates the result obtained from eq. (62). 
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Fig. 9. The radiatively corrected, RG-improved light CP-even Higgs mass is plotted as 
a function of tan /3 for Mg^gy = 1 TeV and m^o = 250 GeV, for two choices of = 
and Xt = 2.4MgugY. See the caption to Fig. 
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In Figs. 1^ and ^, = — 1 TeV, i.e., as Xt = At — /icot P varies, so does At. In 
fact, for m^o ^ m^, the dominant one- loop radiative corrections to m\o depend 
only on Xt and Mg^gy [see eq. (|BUD], so that for fixed Xt, the /i dependence of m/jO 
is quite weak. This is illustrated in Fig. For values of m^o ~ 0{mz), the 
dependence is slightly more pronounced (although less so for values of tan/3 ^ 1) 
as illustrated in Fig. |^. We also display rriho as a function of Mg^gy for a number 
of different parameter choices in Fig. ^. In Fig. |^, we exhibit the tan f3 dependence 
of rrihO for two different choices of Xt. Again, we notice that our approximate 
formula [eq. which is depicted by the dot-dashed line, does remarkably well, 

and never differs from the numerically integrated RG-improved value (solid line) 
by more than 1.5 GeV for Mg^gy < 2 TeV and tan/3 > 1. 

In summary, when the algorithm given by eqs. ( |H^ ) and (|H^) is applied to 
the leading log one-loop corrections plus the leading terms resulting from squark 
mixing, the full (numerically integrated) RG-improved value of rriho is reproduced 
to within an accuracy of about 2 GeV (assuming that supersymmetric particle 
masses lie below 2 TeV). The methods described above also yield accurate results 
for the mass of the heavier CP-even Higgs boson, rriHO . The approximation to the 
radiatively corrected charged Higgs mass is slightly less accurate only because the 
leading rrit enhanced terms are not as dominant as in the neutral Higgs sector.^ 

6 Implications of the Radiatively Corrected Higgs Sector 

Using the results of sections 4 and 5, one can obtain the leading radiative 
corrections to the various Higgs couplings, and proceed to investigate Higgs phe- 
nomenology in detail. Here, I shall describe the procedure used to obtain the 
Higgs couplings and briefly indicate some of the consequences. To obtain radia- 
tively corrected couplings which are accurate in the one-loop leading logarithmic 
approximation, it is sufficient to use the tree-level couplings in which the parame- 
ters are taken to be running parameters evaluated at the electroweak scale. First, 
I remind the reader that tan/3 and m^o are input parameters. Next, we obtain the 
CP-even Higgs mixing angle a by diagonalizing the radiatively corrected CP-even 
Higgs mass matrix. With the angle a in hand one may compute, for example, 
cos(/3 — a) and sin a. These results can be used to obtain the Higgs couplings to 
gauge bosons [eq. (|1^)] and fermions [eq. (|1^)]. Finally, the Higgs self-couplings 
[see Appendix A] are obtained by making use of the Aj evaluated at the electroweak 
scale. The end result is a complete set of Higgs boson decay widths and branching 
ratios that include one-loop leading-log radiative corrections. 

®The approximation to the radiatively corrected charged Higgs mass can be improved by 
including sub-dominant terms not contained in the formulae given in this paper; see Ref. |3l|] for 
further details. 
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The Higgs production cross-section in a two-Higgs-doublet model via the pro- 
cess e^e~ — > Z — i> ZH^{Zh^) is suppressed by a factor cos^(/5 — a) [sin^(/3 — a)] 
as compared to the corresponding cross-sections in the Standard Model. At tree- 
level, we know that the decoupling limit applies when m^o ^ mz- In fact, the 
approach to decoupling is quite rapid as indicated in eq. (pH]). For m^o ^ 2m^, 
the couplings of to vector bosons and to quarks and leptons are phenomenolog- 
ically indistinguishable from those of the Standard Model Higgs boson. Including 
radiative corrections does not alter this basic behavior, although one finds that 



cos2(/? 



a 



more slowly as the radiative corrections become more significant. 



When radiative corrections have been incorporated, new possibilities arise which 
did not exist at tree-level. One example is the possibility of the decay h*^ A^A^, 
which is kinematically forbidden at tree-level but is allowed for some range of 
MSSM parameters ||26| , ^ . We can obtain the complete one-loop leading-log ex- 
pression for the h^A^A^ coupling (assuming m^o ^ mz) by inserting the one-loop 
leading-log formulae for the Aj into eq. (|A.1| ) [0 
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(69) 



If kinematically allowed, /i" — * Al^ Al' would almost certainly be the dominant de- 
cay mode. However, the LEP experimental lower bound on m^o now lies above 
0.5(m/io)max — 62.5 GeV. Thus, the region of parameter space where the decay 
Itl^ Al^ is kinematically allowed is no longer viable. The possibility of mea- 

suring the Al^ Al^ couplings at a future e^e~ linear collider by detecting double 



Higgs production has been discussed in Ref . |^ . Unfortunately, the prospects are 
poor due to low cross-sections and significant backgrounds. 

For the heavier Higgs states, there are many possible final state decay modes. 
The various branching ratios are complicated functions of the MSSM parameter 
space [|^. For example, a plot of the branching ratios of with the leading one- 
loop radiative corrections included, can be found in Ref. This plot indicates 
a rich phenomenology for heavy Higgs searches at future colliders. The precision 
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measurements of Higgs masses and couplings will be one of the primary tasks of 
the LHC and future lepton-lepton colliders |^ , |47[| 

Although the possibility of a light Higgs discovery at LEP still remains, the 
effects of the radiative corrections may be significant enough to push the Higgs bo- 
son above the LEP-2 discovery reach. In this case, the discovery of the Higgs boson 
will be the purview of the LHC. Of course, if low-energy supersymmetry exists, 
then LHC will also uncover direct evidence for the supersymmetric particles. In 
this detailed examination of the Higgs sector, with precision measurements 

of the Higgs masses and couplings, will provide a critical test for the underlying 
supersymmetric structure. Unlocking the secrets of the Higgs bosons will help re- 
veal the mechanism of electroweak symmetry breaking and the nature of the TeV 
scale physics that lies beyond the Standard Model. 
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Appendix 



A Three-Higgs Vertices in the Two-Higgs Doublet Model 



In this Appendix, I list the Feynman rules for the 3-point Higgs interaction in 
the most general (nonsupersymmetric) two-Higgs doublet extension of the Stan- 
dard Model, assuming that the Higgs sector conserves CP. The Feynman rule for 
the ABC vertex is denoted by ig^BC- -^^^ completeness, i?-gauge Feynman rules 
involving the Goldstone bosons [G^ and G") are also listed. 

Interactions involving physical Higgs bosons depend in detail on the parameters 
of the Higgs potential specified in eq. (|T]). 
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2mw r 
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9 
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C/3- 



(Al) 



where I have used the notation 



A3 = A3 + A4 + As . (A.2) 

It is interesting to note that couphngs of the charged Higgs bosons satisfy relations 
analogous to that of mH± given in eq. (|). 

The Feynman rules for three-point Higgs vertices that involve Goldstone bosons 
exhibit much simpler forms 

9h°G'>G° = ^^^ho sin(/3 - a) , 
9h»g°g^ = cos(/3 - a) , 

dh^G+G- = 9h0G0G° ' 

9h°g+g- = 9hogogo ' 
9h'>A0G" = TT^i^hO - m\o) cos(/3 - a) , 

9H0A0G0 = 7r~i^H0 " m\o) sin(/3 - a) , 

9h<^H±G^ = - ^ho) cos{(3 - a) , 

9aoh±g^ = ^^^^H± - ^Ao) ■ (A.3) 



w 



In the rule for the A'^H^G^ vertex, the sign corresponds to entering the vertex 
and leaving the vertex. 

One can easily check that if tree-level MSSM relations are imposed on the A^, 
Higgs masses, and angles a and j3, one recovers the MSSM Feynman rules listed 
in Appendix A of Ref. 
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B Renormalization Group Equations 



In this Appendix, I have collected the one-loop renormalization group equa- 
tions (RGEs) that are needed in the analysis presented in this paper PU| , ^ , P7 
Schematically, the RGEs at one-loop take the form 



where t = \n jjL^ 



(B.i: 



where ^ is the energy scale, and the parameters pi stand for the Higgs boson self- 
couplings Aj (i = 1 . . . 7), the squared Yukawa couplings h?j {f = t, b and r; the two 
lighter generations can be neglected), and the squared gauge couplings (j =3, 2, 
1) corresponding to SU(3) xSU(2) xU(l) respectively. The gj are normalized such 
that they are equal at the grand unification scale. It is also convenient to define 



9 = 92 



9 



IB.21 



where g and g' are normalized in the usual way for low-energy electroweak physics, 
i.e. ta.n9]y = g'/g. 

1 now list the /3-functions required for the analysis presented in this paper. 
Two cases will be given, depending on whether /i is above or below the scale of 
supersymmetry breaking, Mg^gy- 
1. /i > MgusY 
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(B.3) 



Here A^g = 3 is the number of generations, Nh = 2 is the number of scalar doublets, 
and the Higgs-fermion Yukawa couplings are given by 
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(B.5) 

In writing down the RGEs for the Higgs-fermion Yukawa couphngs in eq. ( [B.5| ), 
I have assumed that the Higgs-fermion interaction is the same as in the MSSM; 
namely, $i [$2] couples exclusively to down-type [up-type] fermions. Moreover, in 
deriving the \x < Mg^gY equations, it was assumed that the effective low-energy 
theory at the scale n includes the full two-doublet Higgs sector (but does not 
include the supersymmetric particles, whose masses are of order Mg^gy). 

Finally, I list the RGEs for the Higgs self-couplings of the general two-Higgs 
doublet model (with the Higgs-fermion couplings as specified above). First, I need 
to define the anomalous dimensions of the two Higgs fields: 
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[B.6) 



where the sum over i is taken over three generations of quarks (with Nc = 3) and 
leptons (with Nc = 1). The /5-functions for the Higgs self-couplings in the gen- 
eral CP-conserving non-supersymmetric two-Higgs-doublet model (with the Higgs- 
fermion couplings as specified in section 2) are given by 
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